Screening in semiconductor nanocrystals: Ab initio results and Thomas-Fermi theory 
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A first-principles calculation of the impurity screening in Si and Ge nanocrystals is presented. We 
show that isocoric screening gives results in agreement with both the linear response and the point- 
charge approximations. Based on the present ab initio results, and by comparison with previous 
calculations, we propose a physical real-space interpretation of the several contributions to the 
screening. Combining the Thomas-Fermi theory and simple electrostatics, we show that it is possible 
to construct a model screening function that has the merit of being of simple physical interpretation. 
The main point upon which the model is based is that, up to distances of the order of a bond length 
from the perturbation, the charge response does not depend on the nanocrystal size. We show in a 
very clear way that the link between the screening at the nanoscale and in the bulk is given by the 
surface polarization. A detailed discussion is devoted to the importance of local field effects in the 
screening. Our first-principles calculations and the Thomas- Fermi theory clearly show that in Si and 
Ge nanocrystals, local field effects are dominated by surface polarization, which causes a reduction 
of the screening in going from the bulk down to the nanoscale. Finally, the model screening function 
is compared with recent state-of-the-art ab initio calculations and tested with impurity activation 
energies. 



I. INTRODUCTION 

Screening in semiconductor nanocrystals is a funda- 
mental issue, the importance of which is mostly due to 
the large amount of technological applications inherent 
in the world of nanostructures. While screening in bulk 
semiconductors is a well known and widely investigated 
subject the phenomenon of screening in nanos- 
tructures is still not fully understood. The presence of 
many papers recently published on this subject shows 
that there is a strong interest in this fieldi&^&^ i 10 ' 11 ! 12 
In particular, it has been shown that the response of 
a nanostructure to an external field is bulklike well in- 
side the structurei^ This is an interesting feature due 
to the local nature of bonding in semiconductors. On 
the other hand, such a behavior is quite different from 
the trend shown, in both experimental results and the- 
oretical calculations, by the nanocrystal macroscopic di- 
electric constant) 12 ! 13 ' 14 ! 15 that is shown to increase quite 
slowly with the size upon going from small nanocrystals 
to the bulk. The connection between the microscopic 
local, bulk-like response and the macroscopic dielectric 
properties can be understood only through the study of 
the surface polarization contribution, which strongly in- 
fluences the nanocrystal properties. 

Preliminary results on this subject have been presented 
in Ref. IS. In the present paper, we investigate the na- 
ture of the screening in semiconductor nanocrystals by 
performing a detailed microscopic analysis. In the first 
part of the paper, a first-principles study of shallow impu- 



rity screening in nanocrystals is described. The charge 
density induced by donors and acceptors in Si and Ge 
nanocrystals is discussed, pointing to the connections to 
both the point-charge case and the bulk limit. It is shown 
that, just like in bulk semiconductors^ the response to 
isocoric impurities in nanocrystals is similar to that of 
a point-charge. A description of the local field effects 
on the induced charge is given showing that the present 
real-space analysis gives a way to distinguish between 
surface and bond polarization contributions. In particu- 
lar, a very interesting comparison with the work of Ref. 
[H allows us to give a physical interpretation of the various 
contributions to the screening. 

In the second part of the paper, a model for the screen- 
ing is proposed. The well known electrostatic model, 
based on the image charge method^ is not valid in the 
neighborhood of the impurity, in that wrong boundary 
conditions are given for the induced potential. We pro- 
pose a generalization to semiconductor nanocrystals of 
the Thomas-Fermi model as originally proposed in Ref. 
[H for bulk semiconductors. In order to better appreciate 
the merits and the shortcomings of the Thomas-Fermi 
theory, an alternative derivation of the model from the 
Hohenberg and Kohn theorem^! is described. In the 
present model, both the correct limit in the neighborhood 
of the impurity and the surface polarization contributions 
are taken into account from the beginning as boundary 
conditions of the Poisson equation. Our results fit well 
the screening diagonal contribution^ for bulk silicon, and 
excellent agreement with a recent first-principles calcu- 
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lation of the screening function^ is obtained. Moreover, 
the model gives a fair prediction of the impurity binding 
energies. 

II. AB INITIO RESULTS 

We have performed ab initio calculations using differ- 
ent kinds of impurities in Si and Ge nanocrystals. The 
study has been based on a plane-wave density-functional 
theory (DFT) framework. The calculations have been 
done with the QUANTUM-ESPRESSO package^ us- 
ing the general gradient approximation with ultrasoft 
pseudopotentials. 19 A vacuum space of at least 6 A has 
been left within the supercell, in order to avoid spuri- 
ous interactions between a nanocrystal and its replicas. 
The convergence with respect to the plane-wave basis- 
set cutoff has been treated with care, and Makov-Payne 
corrections 2 ^ have been added for the computation of 
the charged nanocrystal total energies. The nanocrys- 
tals have a nearly spherical shape, they are centered on 
a Si (Ge) atom, and the surface dangling bonds are sat- 
urated with hydrogen atoms. The undoped structures 
have been relaxed, and the optimized geometries used for 
the doped structures. In this way we study the screen- 
ing due to the electronic response to an external pertur- 
bation, where fixed ionic positions are considered. The 
resulting screening does not take into account the con- 
tribution due to the ionic displacement. This latter is 
negligible when covalent semiconductors are considered. 

The electron density induced by a donor species X z in 
the nanocrystal Si;H m is given by 

"ind = n [{Sii-iXR m ) z } - n [ShH m ] , (1) 

where Z is the impurity net charge (atomic units are 
used throughout this work). A similar expression holds 
for acceptors and for Ge nanocrystals. We have focused 
our study on the spherical average of ni n <j, and on its 
integrated density defined as 

QindM = / n ind (x)dx, (2) 
Jo 

where n ind (r) = 47rr 2 ?i ind (r). 

In Fig. [TJ the results for donor and acceptor isocoric 
dopants in SissE^ are reported. The spherically aver- 
aged induced densities fii n( i(r)/Z [Fig. [TJa)] and the inte- 
grated densities Qi n d(r)/Z [Fig. [TJb)] are shown for the 
donors P + , S +2 and acceptors Al~, Mg~ 2 . It is worth 
noticing that the densities induced by P + and S +2 (solid 
and dashed line, respectively) are almost indistinguish- 
able. This suggests that, at least for donor isocoric im- 
purities, the screening is within a linear-response regime. 
This was already argued some years ago for the bulk, 2 
and our analysis shows that this is valid for nanocrys- 
tals too. From Fig. [T| it is seen that, in the case of the 
acceptors AU and Mg~ 2 (circles and stars), apart from 
a narrow region of space near the impurity (inside the 
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FIG. 1: (Color online) Spherically averaged density ni n[ j(r)/Z 
[panel (a)] and integrated density Qi n d{r)/Z [panel (b)] in- 
duced in Si35H36 by different doping species: S +2 (dashed, 
pink line) and P + (solid, black line) donor doping, Al~ (cir- 
cles, red) and Mg -2 (stars, blue) acceptor doping. 



pseudopotential core, the acceptor curves show a nonlin- 
ear peak that is absent in the donor case), the overall 
structure is similar to the donor case. However, while 
in the region of space after the second main peak from 
the impurity the curve is very similar to the donor case, 
the first main peak is lower and pushed away from the 
impurity site (r = 0). This discrepancy between donor 
and acceptor doping is in agreement with some results 
recently published^ and is related to nonlinearity ef- 
fects. The problem of nonlinear screening in bulk semi- 
conductors has been studied with care in the past within 
a Thomas- Fermi formalism. In Refs. 21.22.23.24, it was 
shown that there is an asymmetry in the response when 
positive or negative charged impurities are considered. 
While for donor screening (Z = +1, +2) nonlinear effects 
are negligible, in the case of acceptor screening the re- 
sponse is quite far from linearity. As we shall see in the 
next section, the screening in semiconductors is mostly 
confined within a screening sphere around the impurity. 
The deviation of the radius of such a sphere (screening 
radius) from the linear case furnishes a realistic check of 
the linearity of the response. From the Thomas-Fermi 
theory, an increase of the screening radius of about 0.39 
and 0.68 a.u. (with respect to the linear case) was ob- 
served for Z — — 1 and —2, against a decrease of only 0.16 
and 0.26 a.u. found for Z = 1 and 2, respectivel y 22 ! 24 
These results show the same trend as in Fig. [TJ confirm- 
ing our first hypothesis of a good linearity of the response 
in the donor case. 

In order to show that the features discussed above for 
both the induced and integrated density are not peculiar 
to silicon, we have also studied the case of Ge nanocrys- 
tals. In Fig. [21 our results for the induced [panel (a)] 
and integrated [panel (b)] electron densities are shown for 
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FIG. 2: (Color online) Spherically averaged density fii n a(r)/Z 
[panel (a)] and integrated density Qi n d{r)/Z [panel (b)], in- 
duced Ge35H36 and Gei9iHi48 nanocrystals. The results 
shown are for P + (dashed, blue line) and As + (solid, black 
line) in Ge3sH36, As + (circles, green) in GeigiHi48. In panel 
(b) we also report the linear response to a point-charge im- 
purity in bulk Ge from Ref. (stars, red). 
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FIG. 3: (Color online) Spherical averaged induced density 
nind(r) [panel (a)] and integrated induced density Qi n d(r) 
[panel (b)], calculated for P + impurities in Si nanocrys- 
tals with increasing size: S135H36 (solid-black line), Sis7H7e 
(dashed-red line), Sii47Hioo (circles-blue) and Sii9iHi48 
(stars-green). The full arrows point to the nanocrystal ra- 
dius (see text for definition). 



P + and As + in Ge3 5 H 36 (dashed and solid line, respec- 
tively) and As + in GeigiHug (circles). For comparison, 
in Fig. [UJb) we also report (stars) a previous ab initio 
result obtained using a truly point-charge impurity for 
computing the linear-response screening in bulk Ge.— It 
is evident that the screening near the impurity is much 
more pronounced for the nonisocoric P + impurity than 
for the isocoric As + . The different core of the dopant 
atom with respect to the host atom gives rise, in the first 
case, to the so called central cell corrections, which give a 
nonlinear contribution. On the other hand, the response 
to (isocoric) As + impurity does not induce central cell 
corrections. It should be pointed out that these correc- 
tions have a finite range, being especially pronounced in 
the region of space before the second main peak in the 
QiTid{r)/Z curve. Instead, the electron charge localized 
on the surface is quite independent of the chemical shift. 

Based on Fig. [2j an interesting analysis of the size de- 
pendence of the screening can be deduced as well. Indeed, 
from a comparison between our results for Ge35rl 36 and 
Gei9iHi48, we argue that, upon increasing the nanocrys- 
tal size, the induced charge close to the impurity rapidly 
converges to its bulk value. This can be better inferred 
from panel (b), showingresults at small r very close to 
those calculated in Ref. U for bulk Ge. 

From the results illustrated so far, we can reasonably 
conclude that (i) the isocoric doping is well approximated 
by a point-charge screening, (ii) the screening is in the 
range of the linear response regime, (iii) the induced den- 
sity rapidly converges to the bulk in the region of space 
close to the impurity. 

Further considerations can be driven from panel (b) 
of Fig. O In Ref. 4, it was shown for bulk Ge that 



the induced charge density [stars in Fig H[b)] receives 
contributions from both diagonal and off-diagonal terms 
of the dielectric matrix in reciprocal space. The diago- 
nal contribution to the integrated density consists in a 
monotonic increase of Q m d{r)/Z , followed by constant 
value. Instead, the off-diagonal contribution, related to 
the so called local field effects, gives a wavy structure to 
the integrated induced density. This undulating behav- 
ior has been discussed by several authors for bulk sys- 
tems, and arises from the polarization of the bonds4^i2£ 
Both the diagonal and off-diagonal contributions are re- 
trieved for the nanocrystals, as is clearly seen in Fig. [21 
but an additional effect due to the polarization of the 
surface emerges. It consists of an electron density ac- 
cumulation around the nanocrystal surface due to the 
dielectric nature of the structure, as is well known from 
classical electrostatics. The surface polarization is re- 
trieved from the present ab initio formulation and can 
be seen as a very special kind of local field effect causing 
the annihilation of Q m d {f) /Z in correspondence with the 
nanocrystal boundary [Figs. [BJa) and^a)]. We conclude 
that in a nanocrystal, local field effects (i.e., not diago- 
nal contributions) show up as both the bond and the 
surface polarization. Nonetheless, while the bond polar- 
ization, which is a bulk effect, gives a negligible contri- 
bution to the screening and the optical properties, the 
surface polarization is extremely important. Indeed, it 
is closely related to the depolarization effects that have 
been shown to strongly modify the optical properties of 
low-dimensional structures ! 26 ! 27 ' 28 

The size dependence of the screening can be inferred 
from Fig. [3J where the spherical average of the induced 
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density [panel a] and the integrated density [panel b] are 
shown for a set of P + -doped silicon nanocrystals with in- 
creasing radii. It comes out that up to a distance of a few 
a.u. from the impurity, both the induced and the inte- 
grated densities are almost independent of the nanocrys- 
tal size. In order to be more quantitative, it is necessary 
to define the nanocrystal radius. We do not follow the 
standard route, that is, the radius of sphere whose vol- 
ume equals the product of the volume per atom in the 
bulk times the number of, say, silicon atoms. Indeed, the 
shortcoming of this definition is that the hydrogens passi- 
vating the surface are left out. Moreover, in a screening 
problem the fundamental objects are the electrons, in- 
cluding those participating in the bonds near the surface. 
On the basis of these arguments, we define the nanocrys- 
tal radius through the equation uqAttR^/S = N, where 
uq is the bulk valence electron density and N = Al + m 
is the total number of valence electrons in SiiH m . The 
nanocrystal radii calculated in this way are shown in Fig. 
[3] as arrows. They have slightly higher values than the 
conventional radii. For instance, in the case of Si35H 36 , 
we retrieve 11.3 a.u. against the conventional value of 
10.4 a.u.. But, as it is seen from Fig. [3)Ja), with the 
present definition all the radii lie almost exactly in the 
middle of the surface charge density. As a final obser- 
vation, we stress that despite the undulation in both 
the induced and integrated densities, the total induced 
charge around the impurity is almost independent of the 
nanocrystal radius. Although, at this level, this concept 
may not appear to be well-defined, we shall see in the 
next section how it allows the construction of a Thomas- 
Fermi model for the screening. 



III. THE THOMAS-FERMI MODEL 

As we have seen in the previous section, the density 
induced by a point-charge impurity in semiconductor 
nanocrystals consists of several contributions. The first 
contribution is due to the reciprocal-space diagonal part 
of the dielectric matrix, the second contribution comes 
from the Si-Si bond polarization, and the third is from 
the nanocrystal surface polarization. In this section, we 
illustrate a Thomas-Fermi model, which reproduces with 
great accuracy both the diagonal and surface polarization 
contributions. 



A. Classical electrostatics 

As a starting step, we briefly describe a simple elec- 
trostatic model that can be used as a first approxima- 
tion in studying point-charge screening in semiconductor 
nanocrystals. The basic assumption is that, as a response 
to a point-charge potential v ex t(r) = —Z/r, an electron 
polarization charge — Q is induced around the impurity. 
On the other hand, since the total induced charge must 
integrate to zero, a compensating charge Q, assumed as 



uniformly distributed on the nanocrystal surface, must 
be introduced. The resulting potential energy is there- 
fore due to a total charge Z — Q localized around the 
impurity, and to a charge Q distributed on the surface. 
From electrostatics, the potential energy for an electron 



v c (r) 



Z-Q 



- 1, < r < R 
#, r>R, 



(3) 



where R is the nanocrystal radius, as defined in the pre- 
vious section. 

This electrostatic model is supported by the DFT re- 
sults described in the previous section. Indeed, Fig[3] 
(a) shows that the induced electron density has two 
main contributions concentrated around the impurity 
and across the surface. A basic hypothesis for the model 
is that the total charge Q is independent of the nanocrys- 
tal size. This is supported by the fact that, as shown in 
Fig. [21 the induced density near the impurity rapidly 
converges to the bulk value. As is shown below, the 
use of the bulk value for Q is a good approximation for 
nanocrystals. 

However, this simple model is not accurate in the de- 
scription of the screened potential close to the impurity 
site. Indeed, one can easily see that unphysical poten- 
tial values are predicted in the limit r — > 0. In order to 
correct the model, we follow Ref. [H, in which a Thomas- 
Fermi description for bulk semiconductors is given. The 
basic point is the concept of incomplete screening occur- 
ring in a semiconductor. This means that the charge is 
induced only inside a screening sphere of radius R s . Out- 
side the screening sphere, the system behaves as a classi- 
cal dielectric medium having a static dielectric constant 
s s . In a semiconductor, the total induced charge Q is a 
finite fraction of the whole external charge Z introduced 
with the impurity. It is well known that the ratio Q/Z 
closely depends on the static dielectric constant through 
the relation^ 



Q = 1 _i 

Z e„' 



(4) 



The larger the static dielectric constant is, the closer to 
unit is the ratio Q/Z . The induced potential and the spa- 
tial dielectric function are obtained solving the Thomas- 
Fermi equations with appropriate boundary conditions, 
where the screening radius is derived self-consistently. 3 

The Thomas-Fermi theory establishes the linear rela- 
tion between the induced electron density ni n( j and the 
screened potential w scr = u ex t + ^iud inside the screening 
sphere, 



nindM = —T~ bscr(V) - Sfj] , 
47T 



(5) 



where 8fx (the chemical potential) is a constant to be 
determined with the boundary conditions and q is a mul- 
tiplicative constant related to the nanocrystal average va- 
lence electron density^ Outside the screening sphere, the 
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induced charge is zero and the screened potential matches 
the classical expression v c given in Eq. ([3]). The conti- 
nuity of the induced density, which inside the screening 
sphere is given by Eq. ([5]), at r = R s implies [n.ind(?") = 
if r > R s ] that 



Sfi = v c (R s ). 



B. Derivation of the Thomas-Fermi theory 



(6) 



Before going to the application of the Thomas-Fermi 
theory to a nanocrystal, it is of some interest to rederive 
Eq. ([5]) within the DFT framework. From the Hohcnbcrg 
and Kohn theorem, it is known that the total energy of 
an electron system, written as a functional of the electron 
density, can be decomposed into several contributions, 



E [n] = T s [n] + J [n] + E xc [n] + V [n] 



(7) 



where T s is the noninteracting kinetic energy, V is the 
total one-particle potential energy 



V [n] = J drn(r)v(r) 7 



J is the classical Coulomb energy 

,n(r)n(r') 



J[n] 



dr dr 



(8) 



(9) 



and E xc is the exchange-correlation energy density 
functional] i 29 From the variational principle, with the 
constraint that the integrated density gives the total 
number of electrons N, we have the stationary condition 
for the ground-state density n gs (r), 



SE 
Sn 



(10) 



The variable [i has been shown to be the system chemical 
potentia h 17 ' 29 i 30 In order to study the screening, we write 
down the Hohenberg and Kohn equations, for both the 
unperturbed (labeled with 0) and the perturbed (labeled 
with 1) systems 



Sn(r) 
ST s [n] 



Sn(r) 



+ j)o(r) + ii ff [ji (r)] + iJ !:c [no(r)] = /xq(11) 
+ v 1 (r) +v H [n 1 (r)]+v xc [n 1 (r)} = Mi(12) 



Here, v is the one-electron potential, vh is the Hartree 
potential 



v H [n(r)} 



dr 



, n(v>) 



r' - r 



(13) 



and v xc is the exchange-correlation potential. By sub- 
tracting Eq. (TIT]) from Eq. (jT3J) , and assuming linear 
response, we obtain an expression relating the induced 



density nj n d = n\ — no with the screened impurity po- 
tential w scr . This potential is defined as the sum of the 
external perturbation v cx t = t>i — vq, the induced electro- 
static potential Svh = vh [Sn] , and the induced exchange- 
correlation potential 



SVjcc = / dr 



5n(r) 



nind(r)- 



(14) 



The final result is 



dr' 



S 2 T s [r 



Sn(r)Sn(r') 



<id(r') > +u 8C r(r) =6(1. (15) 



Before showing how and under what conditions Eq. ([5]) 
can be derived from Eq. (|15p . it is worth mentioning 
that this derivation is equivalent to the standard linear 
response theory, in which the perturbation theory is ap- 
plied to the self-consistent Kohn-Sham equations^ 

Equation (Tj"5"|) is an integral equation in which the 
exact form of the kernel is unknown. The Thomas- 
Fermi approximation consists in approximating the ki- 
netic functional to the free-electron gas case. Therefore, 
we have^ 2 - 



S 2 T s [r 



5n(r)6n(r') 



10 

= a—n {r) 



^(r-r'), (16) 



so that the integral condition Eq. (f]~5|) reduces to the 
following algebraic expression: 



«ind(r) = ~ U scr(r)] , 



where 



10a 



n (r) 



1/3 



1/2 



(17) 



(18) 



and a — 3 5 / 3 7r 4 / 3 /10. In the standard Thomas-Fermi 
model no(r) is approximated as the constant spatial av- 
eraged valence electron density, which is the number of 
valence electrons in a unit cell divided by the unit cell 
volume. However, anticipating the results for Si and Ge 
to be discussed below, we have found that this approxi- 
mation is not so important. What turns out to be crucial 
is the choice in Eq. (fTS)) for the second functional deriva- 
tive of the kinetic energy. Indeed, writing it as a function 
of only the difference r — r' corresponds to neglecting the 
local field effects, in that the response is independent of 
the coordinates in which an external point-charge impu- 
rity is located. This is equivalent, in a bulk system, to 
only considering the diagonal component of the recipro- 
cal space dielectric matrix. This argument is about the 
kinetic functional, and does not regard the surface polar- 
ization effects, which are taken into account in the model 
through the use of suitable boundary conditions for the 
Poisson equation. We know that local field effects enter 
the theory through the kinetic functional derivative, as 
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well as through the boundary conditions of the Poisson 
equation. The Thomas-Fermi theory neglects the first, 
but retains the second of such effects. In the bulk, only 
the first effect occurs, due to the polarization of bonds, 
and a link with the diagonal screening of the analysis of 
Ref. |j can be traced. The Thomas-Fermi theory gives 
results in fair agreement with the linear-response results, 
when the off-diagonal elements of the dielectric matrix 
are neglected. 

In order to give a first indication on the validity of 
the Thomas-Fermi approximation in nanostructures, we 
compare the chemical potential Sfi calculated from DFT 
with the Thomas-Fermi prediction. Within DFT, the 
chemical potential of a given structure is calculated as 



C. Model description and results 

The Poisson equation, together with Eq. ([5]), gives the 
Thomas-Fermi equation for an external point-charge Z 
located at the nanocrystal center, 



vV 



AirZS(r) + q 2 [v acx (r) - Sfi] 



It can be written as 



vV 



with the condition 



q 2 [v SCI (r) - A] 



(19) 



lim [rv SCI (r)] = -Z. 



(20) 



(21) 



(22) 



where / = Eg — Es and A = Es — E$ are, respectively, 
the structure ionization potential and the electron affinity 
(Es is the total energy)^ We have calculated the chemi- 
cal potential variation Sfi for Si35H36, for the isocoric P + 
impurity, obtaining Sfi = —2.79 eV from Eq. (119|) and 
Sf.i = —2.62 eV from the Thomas-Fermi Eq. © (for the 
determination of the parameters, see below). Consider- 
ing all the approximations involved, the agreement of less 
than 0.2 eV can be considered very good. 



J 



The leading boundary condition is given by the continu- 
ity of the potential on the screening sphere, 



v scr: (R s ) = v c (R s ). 



(23) 



Solving the Poisson equation, we find that the Thomas- 
Fermi expression for the effective spatial dielectric func- 
tion, defined as the ratio between external and screened 
potential, is 



e(r) = 



{sinh [q(R s r)} + rq} + 
\ ± Z R Z j ' 



< r < R s 
R s <r <R. 



(24) 



Moreover, from the continuity of the electric field on the 
screening sphere, a relationship between the product qR s 
and Q can be obtained, 



sinh(gi? s ) 
qRs 



l-Q/Z 



(25) 



In the following, we shall use Eq. (|25|) for determining 
the variable q for each nanocrystal, once the total density 
Q and the screening radius R s are known. 

The determination of the Thomas-Fermi parameters 
Q and R s is not obvious. We propose to use the bulk 
Q, and calculate the screening radius from the ab initio 
nanocrystal-induced density. As we shall see below, the 
R s dependence on the nanocrystal size is very weak, and, 
a posteriori, we can say that the use of a constant value 
for the screening radius does not significantly change the 
results. For determining Q in the bulk, it is necessary to 
consider Eq. (f2"4"| in the limit R — > oo. It is easy to see 



that 



huik{r) 



{^^{^[q(Rs-r)]+r q }} \ 0<r<R s 



r >R S 
_ (26) 

This is the equation first derived in Ref. y. Since the 
material behaves like a dielectric medium beyond the 
screening radius, we get e s = 1/(1 — Q/Z) from the above 
equation, in agreement with Eq. ([4]). 

From this equation we have Q = 0.91 and 0.93 for Si 
and Gc, respectively, corresponding to the bulk static di- 
electric constants e s — 11.4 and 14.3. The choice for Q 
in the case of Ge is motivated by the fact that it is the 
diagonal contribution to the integrated density, in the 
limit r — > oo^ Also, the corresponding e s is very close 
to the Walter and Cohen pseudopotential calculations of 
dielectric function (they used e s = 14 for bulk Ge).— 
This allows us to show that the present model fits well to 
a standard RPA calculation. Finally, we want to point 
out here that the experimental value e s = 16 is not far 
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TABLE I: Thomas-Fermi parameters calculated for several Si 
and Ge nanocrystals and for bulk Si. The nanocrystal radius 
R, the screening radius R a , and the parameter q, as computed 
from Eq. (|25p . are reported. All the quantities are in atomic 
units. 





R 


Rs 


q 


Si35H36 


11.31 


5.58 


0.836 


Si87H76 


15.17 


5.39 


0.866 


Sii47Hioo 


17.82 


5.35 


0.873 


Sil9lHi48 


19.58 


5.36 


0.871 


Sibulk 


oo 


5.33 


0.875 


( • II- 


11.92 


5.90 


0.84 


Gei9iHi48 


20.62 


5.61 


0.882 



1.2 f 1 1 1 =i r 1 1 1 =i 1.2 




' 10 " 20 10 

r (a.u.) r (a.u.) 



from the value we use. We have checked that small differ- 
ences in the static dielectric constant do not significantly 
change our results. 

It remains now to fix the screening radius R s . This can 
be done starting from the DFT-induced charge illustrated 
in the previous section. One can define R s as the radius 
corresponding to an integrated induced charge equal to 
Q. Namely, the screening radius is given by the solution 
of the equation 

Qind(Rs) = Q. (27) 

However, because of the wavy behavior shown by the in- 
duced charge, there can be several solutions to Eq. (|27| . 
Among the possible solutions, a choice can be done start- 
ing from physical motivations; in particular, we expect 
that the screening radius would take a value close to 
(actually slightly greater than) the impurity host bond 
length. Indeed, the whole displaced charge has to be in- 
cluded in the model, and it is likely that, after giving off 
the wavy contribution, the induced charge extends to the 
nearest-neighbor distance.— 

In Table [I] we report the calculated values for 
the Thomas-Fermi parameters, for several Si and Ge 
nanocrystals. We have also included in this table the val- 
ues for bulk Si that have been obtained with a calculation 
on a supercell comprising 512 atoms. The important re- 
sult coming from Table U is that R s is quite independent 
of the nanocrystal size (apart from the case of S135H36, 
very small variations are seen, being less than 0.1 a.u. 
going from Sis7H7 6 to the bulk Si) . We can say that this 
independence of R s of the system dimension is a mani- 
festation of the local nature of the response as originally 
suggested in Ref. |9|. 

It is important to point out that we use a different 
procedure from that in Ref. [H in the determination of the 
parameters. In that case, q was calculated starting from 
the material Fermi energy, while the screened radius was 
a derived variable, calculated from Eq. ([25]) . At variance 
with that procedure, we directly calculate the screening 
radius from first-principles results and then use Eq. (125|) 
to get q. The values of R s and q obtained in Ref. |_3| were 



FIG. 4: (Color online) Integrated induced density calculated 
for bulk Si [panel (a)] and Ge [panel (b)] using both DFT 
and the Thomas-Fermi model. In panel (a), Thomas- Fermi 
(full, black line) and DFT results (circles, green) for SiigiHi48 
nanocrystal. In panel (b) , Thomas- Fermi (full, black line) and 
linear response results (stars, red) of Ref. [j. The diagonal 
component of screening is also shown (dashed, blue line). 



R s = 4.28 and q — 1.10 for silicon and R s = 4.71 and q = 
1.08 for germanium. Although the present estimations of 
R s (q) are bigger (smaller) than those of Ref. H, we shall 
see below that the dielectric functions for both the bulk 
and nanocrystals are in excellent agreement with DFT 
and empirical pseudopotential results. 

Before presenting the results for the dielectric func- 
tions, let us look at the integrated densities in bulk silicon 
and germanium. The results are shown in Fig. 01 where 
the solid line comes from the Thomas-Fermi model using, 
for silicon, the calculated bulk parameters and for germa- 
nium those corresponding to GeigiH^g (see Table|T|. In 
the case of Si [panel a], we compare the model with our 
nanocrystal DFT results. For Ge [panel b], a comparison 
is done with the linear-response calculation of Ref. H|. We 
note from Fig. [5] that the characteristic wavy structure 
that is typical of the full response (stars) is absent. More 
importantly, the agreement with our Thomas-Fermi cal- 
culation is very good. In other words, the model is per- 
fectly capable of reproducing the response without local 
field effects. This is also true for silicon, as shown in Fig. 
0Ja). Here it is clearly seen that the Thomas- Fermi curve 
lies in the middle of the DFT curve, as in the case of ger- 
manium. We can therefore reasonably conclude that the 
present model gives results consistent with an ab initio 
calculation in which local field effects due to the bond 
polarization are either neglected or play a minor role. 

Within the present model, the wave-vector-dependent 
screening dielectric function of a bulk semiconductor 
takes the form 



q 2 (1 - Q/Z) sin (kR s ) /kR s + k 2 ' 
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0.5 

k (a. u.) 



0.5 

k (a. u.) 



FIG. 5: Wave- vector-dependent screening dielectric function 
for bulk Si [panel (a)] and bulk Ge [panel (b)]. The dashed 
lines are the RPA empirical pseudopotential calculation of 
Ref. |33| . The solid line are the Thomas-Fermi results from 
Eq. (HHJ. 



FIG. 6: (Color online) Screening function e(r) for S135H36. 
Solid (black) line: Thomas-Fermi model with the parameters 
of Table U dotted (blue) line: the same but with the pa- 
rameters of Ref. y; dashed (red) line: real-space ab initio 
calculation of Ref. nil. 



In Fig. [51 the results are shown for both bulk Si and 
Ge. A comparison with the empirical pseudopotential 
calculation of Walter and Cohen^ (dashed line) indicates 
that the model parameters listed in Table U give a very 
good agreement in the bulk limit. 

In order to see how the present model performs in 
the case of nanocrystals, a comparison has been done 
with recent real-space ab initio calculations. In Fig. [6l 
the Si35H36 effective spatial screening function calculated 
with the Thomas-Fermi model and the ab initio results 
taken from Ref. 03 are shown. For the Thomas-Fermi 
model, we have used both the derived values of R s and 
q given in Table IJ (full line) and those of Ref. (dotted 
line). It is seen that in this last case the agreement with 
the ab initio results of Ref. [To| is not perfect, particu- 
larly for the height and the position of the main peak. 
However, considering that there is no need of doing any 
calculation for getting the parameters of Ref. H, depend- 
ing on the application at hand, this may surely be a first 
approximation. In any case, it is interesting to see that 
although the calculations of Ref. U3 include the exchange 
and correlation terms, the agreement between our result 
(full line) and the ab initio one (dashed line) is very im- 
pressive. We want to remark once again that the most 
important part of the screening is due to electron charges 
accumulated both near the surface and close to the im- 
purity. 

We also performed a calculation of the donor impu- 
rity binding energy in order to compare the Thomas- 
Fermi model and DFT results. Within a DFT frame- 
work, the binding energy can be calculated as the dif- 
ference -Edft = Id ~ A u between the doped nanocrystal 
ionization energy and the undoped nanocrystal electron 
affinity i 34 i 35 The binding energy can also be calculated 



TABLE II: Impurity activation energies calculated with both 
the ab initio method (-Edft) and the Thomas-Fermi model 
(E m ) for several silicon nanocrystals. The DFT data are for 
the donor isocoric impurities P + . All the values are in eV. 





-Edft 


E m 


S135H36 


3.05 


2.84 


Sl87H76 


2.27 


2.11 


Sil47Hioo 


1.91 


1.72 


Sii9iHi4g 


1.79 


1.60 



using the perturbation theory as 



En 



e(r)r 



(29) 



where s(r) is the effective spatial dielectric function and 
?p c is the first empty state of the undoped nanocrystal. 
Using the undoped nanocrystal DFT wave function ip c 
and the effective dielectric function calculated with the 
model, we can compute the binding energy and compare 
it with the DFT results. In Table HH the bindin g-encrgy 
results are shown. The interesting thing is that there is 
a systematic difference between the DFT results and the 
values estimated from the model. This size-independent 
contribution is quite small, about 0.2 eV for Si. It can 
be due to both the bond polarization effects and the ex- 
change and correlation contributions, which are not taken 
into account in the present model. 
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IV. CONCLUSIONS 

Screening in covalent semiconductor nanocrystals has 
been studied using both advanced ab initio methods and 
a Thomas-Fermi model. Our DFT calculations have 
shown that isocoric donor dopants essentially behave as 
a point-charge giving an induced charge that agrees well 
with the linear-response approximation. Comparing the 
induced integrated densities with those obtained in Ref. 
0, we have been able to conceptually isolate the features 
directly related to local field effects. It is, at least in prin- 
ciple, not easy to distinguish the several contributions to 
the screening due to local fields. However, in the specific 
cases we have analyzed, there is an indication of the fact 
that surface polarization is the dominant local field effect 
in semiconductor nanocrystals. This was proved combin- 
ing the Thomas-Fermi theory including electrostatics of 
surface polarization with the result that the screening 
function agrees very well with state-of-the-art ab initio 
calculations. 

It is worth mentioning that local field effects related 
to surface polarizations may have dramatic consequences 
also in the optical response, particularly for anisotropic 
structure such as wires and ellipsoidal dots. Indeed, it 
has been shown for both silicon^ 7 , and germanium^ wires 
that the optical frequency-dependent absorption function 
is strongly suppressed for light polarized perpendicularly 
to the wire axis. This suppression does not come out 
when local fields are neglected. From a classical point of 
view, these local fields are again dominated by surface 
polarization, as originally recognized in interpreting the 
polarization anisotropy of porous silicon^ Although the 
present model cannot be generalized to the case of the 
response to an external field, it is of interest to note that 
in both cases, that is, point-charge screening and external 
field, surface polarization plays a fundamental role at the 



nanoscale. 

The first-principles calculations presented in this work 
have shown that the concept of bulklike response to an 
external perturbation introduced in Ref. is valid also 
in a screening problem. Although the induced electron 
density contains some oscillations, we have shown that 
it is possible to define a screening radius that is prac- 
tically independent of the nanocrystal dimension. The 
progressive reduction of the screening action is due to 
the surface polarization. In this paper we have shown 
how this contribution may be described with simple elec- 
trostatics. The important point is that now we have an 
analytical expression for the effective screening function 
whose validity for both the nanoscale and bulk has been 
proved. 

A final note should be devoted to a possible general- 
ization of the model to the case of off-center impurities. 
This could be done, at least in principle, by considering 
the matching of the Thomas-Fermi solution to the full 
image potential of an off-center point-charge. However, 
beyond the technical difficulties in doing that, one should 
consider that the method fails, particularly for very small 
nanocrystals, when the impurity is at a distance from the 
nanocrystal surface that is less than or comparable to the 
screening radius. 
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